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ABSTRACT 


With  the  increasing  application  of  power  electronics  to  shipboard  electric  power 
systems,  there  is  a  need  for  studying  power  systems  in  an  object-oriented,  graphical 
environment  in  order  to  examine  the  interrelations  of  the  various  system  components. 

This  thesis  presents  a  method  of  analy2dng  shipboard  electric  power  systems  in  a 
graphical  environment  The  simulation  method  used  allows  for  real-  time  observation  of 
the  various  power  system  measurements  of  interest 

The  machine  equations,  including  those  for  voltage,  torque,  and  rotor  angle  are 
developed  which  permit  analyzing  the  machine  in  any  chosen  reference  frame  depending 
on  the  study  being  undertaken.  Regulators  for  both  speed  and  field  voltage  are  developed 
and  the  stabilit}^  of  the  speed  regulator  is  examined.  A  means  of  modeling  a  load  for  a 
synchronous  generator  is  then  shown  which  accurately  interrelates  the  loading  of  the 
generator  aiid  the  frequency  and  voltage  output  of  the  machine.  This  load  is  then 
connected  to  the  synchronous  generator  and  two  different  scenarios  are  examined, 
including  a  fifty  percent  load  change  and  a  fault  study.  Graphs  showing  the  output  of 
these  studies  are  presented. 
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L  SYNCHRONOUS  MACHD«:  MODEL  DEVELOPMENT 


The  primary  means  of  shipboard  electric  power  generation  is  with  synchronous 
generators  driven  by  steam  turbines.  The  electrical  and  electromechanical  behavior  of  the 
synchronous  generator  can  be  accurately  predicted  from  the  equations  for  a  three-phase 
salient-pole  synchronous  machine. 

The  rotor  of  the  synchronous  generator  is  equipped  with  a  field  winding  and  one  or 
more  damper  windings.  The  rotor  vrindings  usually  all  have  different  electrical 
characteristics.  Besides  the  differing  electrical  characteristics,  the  salient-pole  machine 
exhibits  a  magnetically  asymmetrical  geometry  due  to  the  varying  air-gap  path  for  the 
magnetic  flux.  This  magnetic  asymmetry  creates  time-varying  self-  and  mutual-inductance 
terms  in  the  voltage  equations  for  the  stator  variables.  By  transforming  the  equations  into 
a  reference  frame  which  rotates  at  die  same  augulai'  velocity  as  die  rulor,  the  magnetic  flux 
path  remains  constant,  thereby  eliminating  the  time-varying  inductance  terms. 

Tlus  chapter  develops  the  voltage  and  electromagnetic  torque  equations  for  the 
synchronous  generator  in  teims  of  machine  variables  and  then  converts  these  equations  to 
the  rotor  reference  frame  via  Park's  transformation  [1].  The  developed  model  assumes  a 
magnetically  linear  system  tliereby  neglecting  the  effects  of  saturation.  Saturation  could 
be  accounted  for  with  the  addition  of  another  simulation  segment  for  the  model,  however 
saturation  data  for  the  machines  studied  was  unavailable. 
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A.  VOLTAGE  EQUATIONS  IN  MACHINE  VARIABLES 

Voltage  equations  for  the  synchronous  generator  will  be  developed  from  a  simple 
two-pole,  three-phase,  wye-connected,  salient-pole  synchronous  machine  shown  in  Figure 
1.  The  number  of  poles  of  the  particular  mtu^hine  can  be  easily  changed  as  this  is  one  of 
the  variables  of  the  lesulting  torque  equations.  The  stator  ^vindings  are  assumed  identical 
with  sinusoidally  distributed  windings,  displaced  120®,  with  Ng  C/quivalent  turns  and 
equivalent  resistance  (subscript  s  denoting  stator  variables).  A  field  winding  and  three 
damper  windings  (fd,  kd,  kql,  and  kq2  respectively)  are  installed  on  the  rotor.  All 
windings  are  assumed  sinusoidal  with  applicable  N  equivalent  turns  and  resistance  r 
(subscript  denoting  particular  circuit).  The  quadrature  axis  (q  axis)  and  direct  axis  (d  axis) 

are  used  with  reference  to  the  magnetic  axis  of  the  rotor  circuits.  Voltages  of  the  rotor 
damper  windings  (v^^,  and  for  monitoring  purposes  only,  the  windings 

themselves  are  physically  shorted. 

Voltage  equations  in  machine  variables  [2]  may  be  expressed  in  matrix  form  as 

(1) 

y,jr  =  r,i,^+p'K^  (2) 

where 


^uici 

■  L. 

K 

^abcj 

4. 

2 


and 


Lj^  —  Lb  cos  20, 


L,= 


-iL,-L,cos2(e,-|) 

-iL,-L,cos2(e,+~) 


-iL,-4cos2(0,-|) 

Sit 

L^  +  4-LbCos2(0,  -  — ) 

-—Lj^-Lg  cos  2(0,  +  Tt) 

2 


-iL,-L,cos2(0,+|) 

-  —  L„  -Ljcos2(0,  +7i) 

2 

2it 

Ii,  +  L^-Lacos2(0,  +  — ) 

(4) 


L.= 


L^;COs6,  L^^cosS, 

L^;COs(e,“)  I^t^rCOsCe,-—)  L^^jSinCe,-^)  L^sinCB,-—) 

L^;COs(e,-f-Y)  ^^2C0S(8,+~)  L^^sin(e,+-Y)  L^sin(6,+-^) 


L,*rfSine, 


(5) 


L.  = 


At*/  +  K 


'mkql 


^kqtkql 

0 


^kqlhi2 

A 


0 

0 


(6) 


In  (1)  and  (2),  the  time  differential  operator  isp.  In  (4)  and  (6),  the  subscript  I  denotes  a 
leakage  inductance  term.  The  subscripts  skql,  skq2,  sfd,  and  skd  denote  mutual 
inductances  between  stator  and  rotor  windings. 

The  magnetizing  inductances  are  defined  as 


i,  =|(i-. -i.)  w 

i-  =|('-. +  i.)  (8) 


where  the  term  incorporating  Lj\  refers  to  the  inductance  term  which  is  independent  of 
rotor  position  and  the  term  with  Lg  is  the  inductance  term  which  is  position  dependent. 
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The  rotor  variableii  are  then  referred  to  the  stator  windings  via  the  appropriate  turns 
ratio  giving  the  foliowing  equations 


Z^cos0,  4a  sin  6,  4a  sin  0, 

K=  4,cos(0, -y)  Z^^cos(G,-y)  ZfcjS!n(0, -y)  4^sin(0, -~)  (10) 

Z<^cos(0, +—)  4,cos(0,+Y)  Z^sjn(0, +-~)  4jSin(e,+Y) 


4«i  4?  ^nq  0  0 

j^/  _  4fl  4»2  +  4«  0  0 

0  0  + 

.  0  0  ^«a  4a  +  4a. 


where  the  princed  variables  (i.e.  )  refer  the  appropriate  rotor  variable  to  the  stator 

windings. 

Tlie  voltage  equations  expressed  in  temis  of  machine  variables  referred  to  the  stator 
windings  are 


v,.J_  pl; 

r;+pL; 


(12) 


Note  that  the  voltage  equations  (12)  are  valid  for  positive  direction  of  stator  cmient 
assumed  out  of  the  stator  terminals.  If  motor  action  is  assumed,  the  sign  preceding 
must  be  changed. 


B.  TORQUE  EQUATION  IN  MACHINE  VARIABLES 

The  field  energy  stored  in  the  coupling  field  of  a  synchronous  machine  can  be 
expressed  as 


1  T 


W,  =  -  Lui)i.b.  ■“  (i.b„y  L;i;^ 


(13) 


where 


Lfr-diag[Lij;  Z^2  (1^) 

The  electromagnetic  torque  is  found  by  taking  the  partial  derivative  of  the  coenergy  with 
respect  to  the  rotor  angular  displacement  0,.  The  magnetic  system  is  assumed  linear, 
therefore  the  coeaergy  equals  the  field  energy,  and  upon  taking  the  partial  derivative,  the 
electrical  torque  is 


where  E  is  the  number  of  poles  of  the  machine.  The  torque  is  positive  for  generator  action 
with  the  positive  direction  of  stator  current  assumed  out  of  the  stator  terminals. 

The  torque  and  rotor  speed  are  related  by 

(16) 
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where  J  is  the  inertia  of  tlie  macliine  in  kilograin  meters’  (kg  m’)  or  Joule  seconds^  (J's^). 
The  input  torque  Tj  is  positive  for  a  torque  applied  to  the  shaft  of  the  synchronous 
machine. 

C.  VOLTAGE  EQUATIONS  IN  THE  ROTOR  REFERENCE-FRAME 

The  machine  voltage  equations  (1)  ind  (2)  are  seen  to  contain  flux-linkage  terms 
which  are  dependent  on  rotor  position  0^,  The  analysis  of  the  machine  behavior  is  greatly 
simplified  by  transforming  these  time-varying  inductances  to  constant  inductance  terms. 
This  can  be  realized  if  the  stator  variables  are  transformed  to  the  rotor  reference  frame. 
R.H.  Park  was  the  first  to  incorporate  a  change  of  variables  in  the  analysis  of  synchronous 
machines  [1]. 

A  change  of  variables  which  formulates  a  transformation  of  the  3-phase  variables  of 
the  stationary  circuit  elements  (thus  the  need  to  refer  the  rotor  variables  to  the  stator)  to 
the  rotor  reference  frame  may  be  expressed  as 

V., = k;u  (17) 

where 

f  = 


A. 


(18) 


7 


*‘aba 


fas 

fit 

Jcs. 


and 

cos0^  cos(6^-“)  cos(6, +-^) 

sin0,  sm(0,  — ^)  sin(0,+'^) 

11  1 

.2  2  2 

®r=|,V(4)4+er(0) 


(19) 


(20) 


(21) 


The  term  K'  denotes  the  transformation  matrix  of  stator  variables  (the  subscript-^)  to  the 
rotor  reference  frame  (the  superscript- r).  It  can  be  shown  that  the  voltage  equations  (1) 
and  (2)  are  transformed  into 


''fldOi  ~  P^qiiOs  (^'^) 

=  r;.-  +  p-K-  (23) 


where 


0 


(24) 
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The  flux-linkages  can  be  expressed  in  the  rotor  reference  frame  by  transforming  (9)  to 
obtain 


rv  1  r  k;l,(k;)~‘  K^L'.lr  sr  i 


(25) 


Expanding  the  voltage  equations  given  by  (22)  and  (23)  gives 


(26) 

v'd,=-rA-<ii,X'^  +  p'k'^ 

(27) 

Vo,=-nioAPK 

(28) 

(29) 

(30) 

II 

(31) 

-,/r  ^  4,^  //r  j_  ^ 

(32) 

and  the  expanded  flux  linkage  terras  become 

~  -  +  A»u(  +  'ijl  ■*■  **,2  ) 

(33) 

+  AnJ  (~4  +  J/rf  ) 

(34) 

o 

«« 

II 

1 

o  * 

(35) 

^  t?  1~  1**91  4ii9  (“*9i  **91  ■•■  **92  ) 

(36) 

X  lu,2~  ^1*92**92  ■*■  ^m9  (~*V  ■*■  **91  ■*■  **92  ) 

(37) 

^  /d“  A/d*/rf  ■•■  (“*(ij  +*/!<  ■*■**,4'  ) 

(38) 

^  ]U“  ^IW*IU  ^md  (~*(ir  ■^*/i/  '^^hl  ) 

(39) 
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Expressing  the  voltage  and  flux  linkage  equations  in  terms  of  reactances  rather  than 
inductances  gives  the  voltage  equations  in  matrix  format  (40),  where  the  expression  for 
the  applied  field  voltage  is  introduced 


e 


ft 

4d 


(41) 


and  the  leakage  and  magnetic  reactances  are  combined  as 


(42) 

(43) 

(44) 

X'^2  =  X'^,  +  X^ 

(45) 

(46) 

(47) 

The  zero  sequence  terms  have  been  omitted  since  balanced  conditions  are  assumed. 

The  model  being  developed  will  be  connected  to  a  non-infinite  bus.  This  implies  that 
the  voltage  and  frequency  do  indeed  vary  with  load,  based  on  the  droop  characteristics  of 
the  particular  regulator.  With  this  in  mind,  it  is  necessary  to  solve  for  currents  based  on  the 
bus  voltages,  therefore  the  c  irrents  are  selected  as  independent  variables.  Rewriting  (40) 
with  the  current  and  current  derivative  terms  separated. 
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5 


11 


(48) 

The  state  equation  format  of  (48)  is  used  for  tlie  simulation,  so  the  current  derivatives 
must  be  solved  for.  As  an  example 

V  =  [All  +  [B]pi 
pi  =  [B]-{v-[A]i} 

and,  with  this  format,  the  current  derivatives  arc  solved  as  shown  in  (50) 


12 


Denoting  the  inverse  matrix  elements  as  [B]'  •‘  =  [b]  = 


6m  6, 


'a  ^^12 


-'21 


L*81 


the  desired  differential  current  equations  are  obtained 


K.  =  '■l.K.  +  Vi  + 


0). 


to. 


to. 


-'66  J 


(51) 


Ki  =  ‘bit"; + '•A + ^  e.  -  ^  ‘'/S  -  Or  ^  Cl 


to. 


to. 


to. 


+  633[t;i,^  -  «] 
+  634 


(52) 


PiA  =  ‘..t  A  +  *  +  Or  ^  ij.  -  O,  ^  -  Or  ^  Cl 


+  f>«K'w  - 

+  644 


to. 


(53) 


Pi.  =  +  r,i^.  -  03r  ~  i.  +  to,  i;;  +  to,  ^  ] 

to,  ''  to,  ^  to,  ^ 


+  625 

+  ^6[i'w  - 


(54) 
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+  r,«.  -  CO,  ^  i;  +  CO,  ^  <2.  +  0),  ^  i;;,] 

CUfc  G)a  OJ* 


+  666[y;i-r*y:2] 


(55) 


/)iw  ==  ^>82^1^*  +  +  “<■  ^1*^2] 

03^  0)5  CU^ 

-h  -  XJ},]  (56) 

+  feBatuiw  ”  ^*5] 

Note  the  similar  form  of  the  q-  and  the  d-sjos  integral  current  equations.  These  will 
greatly  simplify  the  SIMULINK  [3]  model 

D  TORQUE  EQUATION  IN  SUBSTITUTE  VARIABLES 

The  torque  equation  (15),  transformed  into  rotor  reference  frame  variables  is 


T.  =  (|)(f  )[4.,<-C  +  i'fi  +  i'uK  - 


(57) 


and  writing  in  terms  of  reactance 


COtf.  _ 

£0*  £0^ 


X  .  x^ 

=  ^andl^  =-^ 

y.. 


(0^ 


£0. 


(58) 
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gives 


H 


V 


+  »’Ci 


(59) 


The  rotor  speed  is  obtained  from  (16)  as 


po), 


(60) 


£.  ROTOR  ANGLE 

The  rotor  angle  of  a  synchronous  machine  is  a  useful  parameter  to  monitor  the  torque 
production  of  the  machine.  It  relates  the  angulai'  displacement  of  the  rotor  and  stator 
magneto  motive  force  axes.  In  equation  form 

5  =  e,  -  e„ 

.  0(0, (^)  -  (o„{5)]4  +  e,(0)  -  e„(0) 

where  §  is  the  rotor  angle  (in  radians),  and  ^  is  a  dummy  variable  of  integration. 

Equation  (61)  can  be  implemented  by  solving  equation  (61)  in  terras  of  the  available 
parameters.  The  rotor  position,  0^  is  available  by  mtegrating  the  rotor  speed.  However, 
the  voltage  angle  9^  must  be  sohzed  for  ao  follows 

4+e„(0)  (62) 
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and 


Vo.  =  +  ^d.  sin  e, 

=  ^jvl  +  vl  cos[e,  +  tan~^(-^)] 

V 

=  \/vi  +  vl  cos[J  CO,  dt  +  8,(0)  +  tan" ‘(•^)] 


(63) 


therefore,  from  (62) 

e„  =  f<i)„(5)  </5  +  e„{0) 

=  [ (I),  *  +  9,(0)  +  tan^'c-^)  (64) 

J  ' - V - - 


giving 

8  =-■  0,  -e,.  =  tan-‘(-^) 

where  a  positive  rotor  angle  is  assumed  for  generator  action. 


(65) 
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n.  SPEED  AND  VOLTAGE  REGULATORS 


For  the  machine  model  to  accurately  simulate  a  shipboard  power  system  it  is 
necessary  to  include  controllers  for  both  the  field  voltage  and  the  prime  mover  torque. 
The  two  principal  controllers  of  a  turbine  generator  set  are  the  automatic  voltage  regulator 
and  the  speed  governor.  The  automatic  voltage  regulator  model  consists  of  voltage 
sensing  equipment,  comparators  and  amplifiers  controlling  a  synchronous  machine  which 
can  be  generating  or  motoring.  The  speed  governor  may  be  considered  to  have  similar 
equipment  but  in  addition  it  is  necessary  to  take  the  turbine  into  account 

A.  SPEED  REGULATOR 
1.  Model  Development 

The  block  diagram  of  a  steam  turbine  and  associated  governor  for  a  typical 
naval  ship's  service  turbine  generator  (SSTG)  is  shown  in  Figure  2  [4]. 

The  reference  speed  input  is  the  no-load  frequency  selected  by  the  operator  of 
the  electric  plant.  The  gain  K  is  used  to  adjust  the  speed  droop  of  the  regulator,  is  the 
time  constant  associated  with  the  servo  and  steam  valve  system  andXp  is  the  time  constant 
associated  with  the  conversion  of  steam  to  torque  ).  The  following  values  were  used 

for  this  study 

x,  0.1  sec 

Xp  0.06  sec 

K  150 

From  this  block  diagram  (Figure  2),  a  SIMULINK  model  was  developed  and  is  shown  in 
Figure  3. 
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As  can  be  seen,  the  speed  regulator  is  designed  to  compare  the  reference  speed  witli  the 
speed  being  produced  by  the  synchronous  mxhine.  This  speed  difference  determines  the 
prime  mover  torque  being  applied  to  the  generator.  When  the  speed  difference  is  no 
longer  changing,  the  prime  mover  torque  is  equal  to  the  electrical  torque  developed  by  the 
generator  (load  dependent).  The  final  machine  speed  (and  therefore  the  machine  output 
frequency)  will  be  at  some  value  below  the  no-load  set-point  frequency.  This  can  be 
demonstrated  by  lefening  to  the  'house'  diagram  shown  in  Figure  4. 

As  the  output  power  of  the  machine  is  increased,  the  operating  frequency  of  the 
machine  decreases  (based  on  the  slope  or  droop  of  the  particular  regulator)  the  operating 
frequency  can  be  restored  to  a  desired  frequency  (i.e.,  60  Hz)  by  the  operator  adjusting  the 
no-load  setpoint  of  the  regulator. 

2.  Regulator  Stability 

As  was  mentioned  in  the  previous  section,  tlie  speed  droop  characteristics  of  the 
regulator  are  determined  by  the  gain  K  of  the  regulator.  The  percent  speed  regulation  is 
defined  by 


%  Speed  Regulation  =  ^ ^  100%  (62) 

“/) 

As  the  gain  K  is  raised,  the  percent  speed  regulation  approaches  100%,  at  wliich 
point  the  regulator  and  therefore  the  regulated  plant  will  become  unstable.  This  stability 
point  was  determined  by  examining  the  root  locus  of  the  model  of  Figure  3.  Ihe 
synchronous  machine  model  was  obtained  by  linearizing  the  machines  current  block 
around  its  stable  operating  point  for  rated  conditions  by  using  the  MATLAB  [5]  fimction 
linnwcL  The  linearized  plant  model  returned  by  linnwd  is  in  linearized  state  space  format 
which  is  then  i  iianged  to  transfer  function  format  with  the  function  ss2tf.  The  transfer 
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Figure  4.  house  Diagram. 
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function  format  of  the  plant,  speed  block,  and  speed  regulator  is  then  reduced  to  a  single 
open  loop  to  study  the  roots  of  the  system  as  the  regulator  gain  K  is  varied.  This 
complete  process  is  demonstrated  in  Figures  5,  5a,  5b,  and  5c.  The  resulting  root  locus 
plot  is  shown  in  Figure  6.  The  operating  constrabt  was  determined  for  the  first  root 
meeting  the  imaginary  axis  for  a  given  K.  The  limit  of  K  was  determined  to  be  153. 
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Figure  Sa.  Current  Block. 


Figure  5c.  Transfer  Function  Development,  Speed  Regulator. 
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B.  VOLTAGE  REGULATOR 


The  voltage  regulator  or  excitation  system  used  in  this  study  is  shown  in  block 
diagram  form  in  Figure  7  [6],  The  parameters  are  defined  as  follows 

and  Regulator  amplifier  gain  and  time  constant 

and  Kg  Saturating  exciter  gain  and  time  constant 

,  and  Kp  Feedback  loop  gain  and  time  constants 
The  SIMULINK  model  developed  from  the  block  diagram  is  shown  in  Figure  8.  The 
regulator  compares  the  bus  voltage  with  the  specified  reference  voltage.  As  can  be  seen  in 
the  Simulink  model,  this  comparison  is  done  in  the  qdO  reference  frame. 

The  reference  voltage  is  computed  as  follows 

^rtfir,nc.=^|vJ'+v/  (63) 

This  equates  to  the  desired  line-to-neutral  voltage  (peak)  as  shown  below 

(M) 


1 _ 

_2 
~  3 

cos8,.  cos(0,  — “) 

sin9,  sin(0^  — — ) 

cos(0,  +~) 

sin(0,+y) 

■  1 

1 

1 

1 

1 

1 

LvJ 

2 

2 

2 

(65) 
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where 


Va=VfaCOS(e„) 

Vi=Vi,COS(e^-y) 

27C, 

V.=VtoCOS(0„+Y) 


(66) 


and 


0,=|a,®4</+e,(O) 

9„=£a),®y+8„(0) 


and  solving  for  Vq  and  v^j 


(67) 


\  =  3V. 


cos0,cos0„  +cos(0,  — ^)cos(0^  — ~)+cos(6r  +-~)cos(0,,  +^)  (68) 


2 

’3^ 


Vd  =“Vt, 


sin0,cos0^+sin(0,  — ^)cos(0,^  — ~)+sin(0,+^)cos(0^  +—)  (69) 


and  using  the  trigonometric  identities 


cos  xcosy+ cos  (.x  — ^)cos  ( y  -”) + cos  (j:+-^)  cos  (>>  +— ) cos  (x  -  y)  (70) 

.  ,  2k.  ,  2n..,2K.  .  2n.  3  .  .  . 

sm  xcosy+sm  {x — — )co3()' — ^^)+sm  (jc+— )cos(y+— )=-sin  (jr-)*)  (71) 

j  3  3  3  ^ 
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gives 


’■,=|vi,|cos(9,-0„) 

(72) 

= v^,  cos[(ai,  -CD,)/ +e,(0)-e„  (0)] 

v<,~|vi,|sin(e,-e.,) 

(73) 

=v^sm[(cD,  -c),)/+e,(0)-e,,(0)] 

At  synchronous  speed  O),  =co,,  and  from  (61)  8=6, -6„,  therefore 

v,  =  Vj,cos(8) 

(74) 

Vj=V4,sin(5) 

(75) 

and 

=  Vvh[cos^(8)  +  sin^(8)]  = 

(76) 

I 
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m.  MACHINE  LOAD 


To  accurately  study  the  machine's  perfoimance  it  is  necessary  to  be  able  to  load  the 
given  machine  in  a  realistic  manner.  Ihe  inductive  nature  of  transmission  lines  and 
induction  motors  dictate  a  need  for  a  resistive-inductive  (RL)  loading  of  the  machine. 

This  section  develops  the  RL-load  used  for  the  synchronous  generator  and  then  the 
SIMULINK  implementation  of  this  loading. 

A.  PASSIVE  RL-LOAD 

The  load  used  for  this  study  was  comprised  of  a  passive  resistive-inductive  load. 
Tuis  loading  is  typical  of  ac  motors,  transmission  lines,  elechic  heating,  and  lighting.  The 
load  was  modeled  iu  a&c-variables  in  order  to  eliminate  the  need  for  additional 
transformations  iabc  qdO).  The  simulation  uses  voltage  as  the  input  and  currents  as  the 
outpvits.  This  requires  the  phase  voltages  to  be  made  up  from  the  network  equations.  The 
loads  also  use  voltage  a^  the  inputs  and  current  as  the  outputs.  Based  on  this  topology  it 
is  advantageous  to  assume  a  Y-connected  set  of  large  resistors  are  connected  across  the 
terminals  of  the  synchronous  generator  [7]  and  the  phase  voltages  at  the  terminals  of  the 
generator  are  computed  from 


(85) 

OO 

N—/ 

Vc=--«(‘cC+U) 

(87) 

where  and  are  the  phase  currents  flowing  into  the  generator  from  the  load  and 

1^,  iti,  and  art  the  phase  currents  flowing  into  tlie  load  from  the  generator.  These 
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currents  are  computed  first  as  and  Jq  and  then  transformed  to  i„,  and  4.  This 
is  repre-sented  in  Figure  9. 

This  load  modeling  is  quite  accurate  as  the  resistance  value  R  is  made  large  (20-50 
pu).  As  the  resistance  is  made  large,  the  current  through  this  branch  becomes  negligible. 

B.  SEVfULINK  LOAD  MODEL 

The  passive  RL-\oad  was  modeled  in  SIMULINK  by  starting  with  the  voltage 
equations  (only  developed  for  phase  a,  others  similar)  from  Figure  10. 

The  difierential  equations  can  be  implemented  as  shown  in  Figure  1 1. 

As  can  be  seen,  the  terminal  voltage  is  dependent  on  the  output  current  .  This 
voltage  dependency  is  the  difference  between  a  finite  and  an  infinite  bus  model  of  electiic 
machines.  An  infinite  bus  analysis  assumes  no  dependency  between  and  4. 

An  attempt  was  made  to  make  the  load  purely  resistive.  SIMULINK  is  unable  to 
deal  with  this  type  of  load  however  since  this  would  incorporate  direct  feed-through 
between  input  and  output  of  the  load.  This  forms  an  algebraic  loop  and  cannot  be 
analyzed.  The  inductance  of  the  transmission  line  overcomes  this  limitation.  The  purely 
resisitve  load,  connected  by  a  transmission  line  (with  inherent  inductance)  to  the  generator 
bus  adds  sufficient  inductance  to  warrant  modeling.  This  inductance  then  incorporates  a 
Ume-delay  between  input  and  output  due  to  the  integral  form  of  the  load  equations. 


(88) 

(89) 

(90) 


Figure  9.  Generator  Terminal  Voltage  Model. 
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Figure  11.  Differential  Equation  Implementation  of  Voltage  Equations. 
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IV.  MACHINE  MODEL  PERFORMANCE 


The  synchronous  machines  performance  is  examined  in  this  section.  The  machines 
parameters  are  given  and  the  machines  response  to  both  a  load  change  and  a  fault 
condition  are  studied. 

A.  MACHINE  PARAMETERS 

The  machine  being  simulated  is  a  4-pole,  60  Hz,  440  V  (line-to-line),  500  kW,  0.85 
power  factor  synchronous  generator.  The  per  unit  parameters  on  the  machine  base  (Z  i^^se 
=  0.1291  O)  are  as  follows 

stator  resistance,  =0.00031 

rotor  field  resistance,  =  0.00226 

q-axis  damper  winding  resistance,  -  0.0224 

d-axis  damper  winding  resistance,  -  0.00673 

stator  leakage  reactance,  Xi^  =  0.0632 

d-axis  magnetizing  reactance,  =  1.4648 

q-axis  magnetizing  reactance,  =  U.7668 

field  leakage  reactance,  X/yy  =  0.1545 

q-axis  damper  winding  leakage  reactance,  X/j^  =  0.0771 

d-axis  damper  winding  leakage  reactai.ee,  X/^y  =  0.3676 

rotor  inertia,  J  =  14.526 

These  machine  parameters  are  taken  from  the  Machinery  Systems  Laboratory  (MSL)  test 
facility  for  a  typical  shipboard  synchronous  generator  [8].  The  MSL  test  facility  is  within 
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the  Machinery  Systems  Engineering  Branch  at  David  Taylor  Naval  Ship  Research  and 
Development  Center. 

B.  SIMULATION:  LOAD  CHANGE 

The  machine  was  initially  loaded  at  an  approximate  power  factor  of  0.8,  with  the  load 
calculated  as  follows 

Machine  data:  440-V  1.5  MVA .  From  Figure  12 


/I  =  =  1.969 

•  '  3  254 

(91) 

7  s  1.969  kA 

(92) 

=129e'“”’mn 

1.969e--'“'"  U 

(93) 

ri=103.2ma  X^- 11.42  mil 

(94) 

With  the  machine  loaded  as  above,  the  expected  steady-state  power  delivered  to  the 
load  will  be  approximately  1.2  MW.  'Fhis  does  not  however  account  for  the  losses  in  the 
stator  windings.  The  machine  tlien  underwent  a  50  %  step  load  change  by  doubling  the 
values  of  r^and  .  This  was  done  in  SIMULINK  by  the  load-change  block  shown  in 
Figure  13.  The  plots  of  the  response  are  shown  in  Figures  14  through  21.  The  data 
obtained  from  the  simulation  was  as  follows: 

FULLY  LOADED  50%  LOADED 

MODEL  MODEL 


Load  power,  MW 

1.1754 

0.5186 

Machine  power,  MW 

1.1792 

0.5151 

Frequency,  Hz 

59.99 

64.32 

Power  Factor 

0.7992 

0.7832 

40 


Figure  13.  SIMULINK  Load  Change  Block. 
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volts  volts 


volts  volts 


?1!QA 


va,  50%  load  decrease 


Figure  16.  50  Percent  Load  Decrease,  va. 
45 


degrees  degrees 


rad/sec  rad/sec 


wr,  50%  load  decrease 


time  (sec) 


Figure  21.  SO  Percent  Load  Decrease,  rotor  speed  (dif). 
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C.  SIMULATION:  THREE-PHASE  FAULT 

The  stability  of  a  synchronous  machine  throughout  a  power  system  following  a  fault 
is  of  major  concern.  A  3-phase  fault  or  short-circuit  rarely  occurs  and  a  3-phase  fault  at 
the  machine  terminals  is  even  more  uncommon;  however,  it  is  stih  instructive  to 
demonstrate  the  ability  to  insert  a  fault  for  a  chosen  time  interval  and  study  the  machines 
dynamic  performance. 

The  machine  was  initially  loaded  at  a  power  factor  of  0.8  as  in  part  B.  A  3-phase 
fault  was  inserted  and  was  left  in  place  for  7  cycles.  The  fault  impedance  (representative 
of  the  bus  work)  is  0.001291  4  j  0.1291 If  the  fault  had  been  allowed  to  remain  on  the 
system  slightly  longer  the  machine  would  have  become  unstable  in  terms  of  not  having 
returned  to  synchronoas  speed  after  removal  of  the  fault  or  the  machine  would  have 
slipped  poles  prior  to  returning  to  synchronous  speed.  It  was  determined  that  a  seven- 
cycle  fault  is  representative  of  a  typical  breaker  reclosure  timing  cycle  which  the  machine 
could  have  been  subjected  to. 

The  three-phase  fault  was  performed  in  SIMULINK  by  the  fault  block  shown  in 
Figure  22.  The  fault  implemented  by  using  a  step  function  block  which  stepped  in  the 
fault  impedance  at  a  particular  time  and  then  another  step  function  block  stepped  in  the 
normal  load  at  the  end  of  the  desired  fault  time.  The  plots  of  the  response  are  shown  in 
Figures  23  through  30. 
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‘':ONCLUSlONS 


A.  ADVANTAGES 

llie  advantages  studying  naval  power  sysiem  components  in  a  graphical 
environment  are  many.  To  be  abl’  to  see,  in  teal  time,  any  Une  or  load  parameter  (i.e. 
current  or  voltage)  makes  the  study  of  system  component  interaction  very  straightforward. 

The  concept  of  "closing  the  loop"  to  make  the  input  and  output  of  the  machine 
dependent  on  one  another,  as  opposed  to  the  iniSnite  bus  concept,  is  easily  implemented 
with  the  graphical  tools  available.  Navy  power  systems  definitely  require  this  non-infinite 
bus  modeling  due  to  their  capacity  difference  as  compared  with  commercial  power  grids. 

B.  ANALYSIS  OF  SIMULINK 

Matlab's  SIMULINK  is  an  excellent  tool  for  the  type  of  studies  conducted  in  this 
thesis.  Model  building  and  integration  into  the  "power  system"  was  easily  done. 

One  of  the  difficulties  with  SIMULINK  included  the  algebraic  loop  problem 
necessitating  an  integrator  in  the  path  to  provide  a  time  delay  element  In  this  particular 
study  this  was  not  a  problem  due  to  the  resulting  differential  equations,  however,  a 
possible  purely  resistive  load  would  be  troublesome. 

On  the  issue  of  speed  of  the  simulations,  naturally  one  is  always  looking  for  a  faster 
simulation.  I  began  my  research  on  a  486-33MHz  machine,  then  graduated  to  a  SPARC-2 
machine  and  realized  a  speed  difference  of  about  two  times.  1  then  moved  to  a  SPARC-10 
and  achieved  a  speed  increase  of  approximately  four  times  over  Ute  SPARC-2.  Typical 
tiraes  to  mn  a  simulation  of  approximately  5-seconds  was  about  4-hours  on  the  SPARC- 
10.  In  comparing  this  to  a  non-graphical  simulation  tool  such  as  Advanced  Continuous 
Simulation  Language  (ACSL)  [9],  a  comparable  model  runs  approximately  in  real  lime  (60 
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seconds  of  nin  time  data  in  approximately  60  seconds).  The  disadvantage  of  ACSL  is,  of 
course,  the  inability  of  the  user  to  observe  system  output  in  real  time  and  the  non-modular 
modeling  of  systems. 

C.  FUTURE  WORK 

The  modeling  of  additional  load  components  including  induction  machines,  another 
synchronous  generator  sharing  the  load,  parallel  RL-loading,  and  adjustable  speed  drives 
could  all  be  performed  in  SIMULINK.  As  the  number  of  integrators  in  the  complete 
model  increases,  however,  the  run  dme  will  also  increase.  In  addition,  the  existing  RL- 
load  can  be  examined  under  such  conditions  as  critical  clearing  time,  first  swing  transient 
stability  limit,  and  various  fault  conditions  such  as  single  line  to  ground  faults,  line-to-line 
faults,  etc. 

In  addidon,  voltage  regulator  stability  can  be  studied  to  find  gain  limits  for  this 
particular  regulator.  Different  regulator  designs  could  also  be  examined. 

Mathworks  has  also  advertised  a  software  package  which  will  take  existing 
SlMUlrlNK  models  and  compile  them;  this  will  supposedly  increase  simulation  speed 
considerably.  This  product  sounds  very  encouraging. 
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